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Abstract 

According to empirical observations, some pattern formation phenom- 
ena in driven many-particle systems are more pronounced in the presence 
of a certain noise level. We investigate this phenomenon of fluctuation- 
driven ordering with a cellular automaton model of interactive motion 
in space and find an optimal noise strength, while order breaks down at 
high(er) fluctuation levels. Additionally, we discuss the phenomenon of 
noise- and drift-induced self-organization in systems that would show dis- 
order in the absence of fluctuations. In the future, related studies may have 
applications to the control of many-particle systems such as the efficient 
separation of particles. The rather general formulation of our model in 
the spirit of game theory may allow to shed some light on several different 
kinds of noise-induced ordering phenomena observed in physical, chemi- 
cal, biological, and socio-economic systems (e.g., attractive and repulsive 
agglomeration, or segregation). 



Noise-related phenomena can be quite surprising and exciting. Therefore, 
they have attracted the interest of many researchers. For example, we mention 
stochastic resonance [1], structural (in)stability [2, 3], noise-driven motion ('Brow- 
nian motors') [4, 5, 6], and "freezing by heating" [7]. The approach proposed in 
this contribution differs from these phenomena. Moreover, since both, the initial 
conditions and the interaction strengths in our model are assumed to be indepen- 
dent of the position in space, the fluctuation-induced self-organization discussed 
later on should be distinguished from so-called "noise-induced transitions" in sys- 
tems with multiplicative noise as well, where we have a space-dependent diffusion 
coefficient which can induce a transition [8] . Although our model is related with 
diffusive processes, it is also different from reaction-diffusion systems that can 
show fluctuation-induced self-organization phenomena known as Turing patterns 
[9, 10, 11, 12, 13, 14]. 

The main goals of this contribution are (i) to qualitatively understand the 
observed noise-induced ordering phenomena in certain self-driven many-particle 
systems and (ii) to derive and investigate a simplified mathematical model for 
them. Some of the properties of this model appear reminiscent of driven (e.g. 
shared) heterogeneous granular systems. For example, one finds stratification 
phenomena such as the segregation of different kinds of grains into layers [15, 16]. 
This mechanism explains, for example, some of the geological formations and ore 
concentrations in the earth. As it may be also used to separate different materials, 
it is interesting to ask for the most efficient separation technique including the 
optimal noise level. A similar segregation phenomenon of lane formation has 
been found in pedestrian counter-streams [17]. This is based on a reduction of 
the interaction strength and related with an increase in the efficiency or average 
speed of motion [18]. 

In the following, we will explore a cellular automaton model which extends 
the one discussed in Ref. [18] by the effects of fluctuations, drift, and asymmetric 
interactions. It may, for example, reflect the interactive one- dimensional motion 
of some driven many-particle systems in a spatial direction perpendicular to the 
main flow direction(s) and to the boundaries. For this purpose, let us imagine the 
example of pedestrian streams in a corridor with two opposite flow directions a. 
We subdivide the one-dimensional space into cells % comparable to the shoulder 
width (diameter) Ax of the pedestrians (particles, entities). Speaking in more 
general terms, we have A subpopulations a with N a entities a distributed over 
I cells % G {1, . . . ,/}. We denote the number of entities in cell % belonging to 
subpopulation a by n". Moreover, we represent the kind of interaction and the 
interaction strength between two entities of subpopulations a and b by a constant 
parameter value P a b- (Generalizations are easily possible.) Finally, we update 
the locations of the entities a according to the following rules: 1st step: Select 
the entity a randomly. If a is located in cell % and belongs to subpopulation a, 
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determine 

S a (j,t)=J2Pa b n b 3 (t)+^(t) (1) 

b 

for j — i and the nearest neighbors j — i± 1, where are random fluctuations 
uniformly distributed in the interval [— p a ,Pa], so that p a denotes the fluctuation 
strength. 2nd step: Move to the neighboring cell % ± 1 with probability 

P a (i± l\i;t) oc max{0,S a (i ± l,t) - S a (i,t)} . (2) 

3rd step: Repeat steps 1 and 2 until the locations of N = J2 a N a entities were 
updated. J^th step: With probability V^°, move all entities a of subpopulation a 
by one cell into the same direction. 5th step: Return to step 1. 

Formula (1) calculates the expected effect of interactions with other entities. 
S a (i,t) is a potential function, see Eq. (4). In the language of game theory, it 
can be called the expected success, since, according to the proportional imitation 
rule (2), an entity a moves to a neighboring cell % ± 1 only if it can increase the 
value of S a . The values P a f, may be interpreted as payoffs in interactions between 
two entities of subpopulations a and b. P a b is positive for attrative, cooperative, 
or profitable interactions, while it is negative for repulsive, competitive, or loss- 
making interactions. The fourth step reflects a bias in the motion of the particles 
of subpopulation a, i.e. a drift velocity. 

We have carried out various simulations with random initial and periodic 
boundary conditions, in order to have a translation-invariant system. Note that 
self-organized pattern formation in such a system implies spontaneous symmetry- 
breaking and a pronounced history-dependence of the resulting state. The typical 
solutions are dependent on the specified payoffs P a b, fluctuation strengths p a , and 
drift velocities V^°. Replacing the asynchronous (random sequential) update of 
steps 1 to 3 by a parallel update yields similar (but less random, i.e. "nicer 
looking") results. Replacing the parallel update of the velocities in step 4 by 
an asynchronous one induces such a high noise level that the system is often 
disordered. 

To obtain a theoretical understanding of our simulation results, we have de- 
rived mean value equations for the densities p a (x,t) = n^(t)/Ax with x = iAx. 
For this purpose, we have derived a master equation and determined the drift- 
und fluctuation-coefficients in the usual way. By second order Taylor expansion, 
the resulting equations can be then approximately cast into the form of Fokker- 
Planck equations [19] 

^ + £ Mx , mx , t)] = D *e^ (3 ) 

where the diffusion coefficients D a increase with the fluctuation strength p a 
in a roughly proportional way. Moreover, D a vanishes when p a vanishes and 
\dS a (x,t)/dx\ is small (as for our homogeneous initial conditions). The exact 
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relation for D a is rather complicated, but not of interest, here. Finally, the drift 
coefficients are given by 



V a (x,t) = V a °+ a ^ 1 with S a (x,t)=Y,Pa bPb (x,t). (4) 

b 

Accordingly, we are confronted with non-linearly coupled Burgers equations, 
which may show a diffusion instability. In the following, we will check whether the 
above mean value equations yield qualitatively meaningful results, i.e. whether 
correlations can be neglected. For this purpose, we will carry out a linear sta- 
bility analysis and compare the theoretical phase diagram with the numerically 
determined one. In the case of two subpopulations a G {1, 2} and V® = V 2 = 0, 
the homogeneous solution p° a = n-/Ai with n\ = N a /I should be unstable with 
respect to fluctuations, if 

p°P n + p° 2 P 22 > D 1 + D 2 (5) 



or 

P1PIP12P21 > (PiPu ~ D 1 )(p° 2 P 22 - D 2 ) . (6) 

Let us first discuss the the symmetric case with p\ = p\ = p, P n = P 22 = P, 
P\ 2 = P 2 i = Q, and vanishing diffusion Di = D 2 = 0. Then, condition (5) reduces 
to 2pP > 0, and condition (6) becomes Q 2 > P 2 . We can distinguish the following 
solutions (see Fig. la, for representatives see Fig. 4 in Ref. [18]): A) If P < and 
q2 < p2 |--_ e _ p < q an( j p < q < _p)] ) a homogeneous distribution p a (x, t) = p° a 

over all sites in both subpopulations is stable with respect to small perturbations 
(which corresponds to disorder). B) If P < (self-repulsion) and Q < (repulsion 
between the subpopulations), but Q < P, we should find segregation (with a 
tendency that all sites are equally occupied, but either by one subpopulation or 
by the other). C) If Q < (repulsion), but P > (self-attraction), we expect 
repulsive agglomeration (i.e. both subpopulations should cluster at different sites, 
with empty sites in between). D) If Q > (attraction) and Q > —P, we should 
have attractive agglomeration (clustering of both subpopulations at the same 
sites, with empty sites in between). Consequently, on the line Q = (P — l)/2 
(i.e. for P = 2Q + 1), we should cross the phase boundary between disorder and 
segregation at P = — 1, the one between segregation and repulsive agglomeration 
at P = 0, and the one between repulsive and attractive agglomeration at P = +1. 
This is, in fact, confirmed by our numerical simulations (see Fig. lb), so that we 
can trust the instability analysis based on the mean value equations. The reason 
for this is the local nature of the interactions. To characterize the different states, 
we have used order parameters of the form 

1 1 _ 1 1 

e (y) = yEfe-^) 2 with Vi = -fY.yi ( 7 ) 

1 i=l 1 1=1 
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to measure the variances of (i) yi = {n\/n\ + nf/nf) (i.e. the deviation from 
a homogeneous occupation of all sites), or (ii) = (n}/nj — nf/nf) (i.e. the 
difference in the degree of occupation by different subpopulations). G(n x + n 2 ) is 
sensitive to (attractive or repulsive) agglomeration (i.e. to clustering with empty 
sites in between), and 9(n x + n 1 ) recognizes, when the two subpopulations tend 
to use different sites (as for segregation or repulsive agglomeration). 
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Figure 1: (a) Theoretical phase diagram of the four qualitatively different patterns 
resulting for the symmetric model without diffusion: A = disordered phase, B 
= segregation, C = repulsive agglomeration, D = attractive agglomeration, (b) 
Order parameters along the line Q = (P — 1)/2 (see solid line in Fig. la), averaged 
over 20 runs after a time period of 20 000(^1+^) update steps with I = 200 cells 
and Ni = N2 = 2000 entities in each of two subpopulations. The theoretically 
predicted phase transitions at P = — 1, P = 0, and P = 1 (see black triangles in 
Fig. la) are clearly visible. 



Let us now focus on the general case with arbitrary payoffs and diffusion. 
While Di,D 2 > increases the threshold for pattern formation in Eq. (5), in 
Eq. (6) it can reduce the threshold for moderate diffusion, while the threshold 
will be higher for large diffusion. We are not surprised that sufficiently large 
diffusion will always give rise to disorder and suppress pattern formation. It 
is interesting though that a medium level of diffusion may cause pattern for- 
mation where the system would otherwise be disordered. Let us focus on the 
example with Pn = —2, P i2 = 2, P 2 i = —2, and P22 = 1, where subpopu- 
lation 2 is repelled from subpopulation 1 and where the self-interaction within 
subpopulation 1 is repulsive as well, while the other interactions are attractive. 
According to conditions (5) and (6), we expect disorder for small fluctuation 
strengths pi, p2- While increasing values of P2 should be counterproductive, in- 
creasing pi should be able to produce pattern formation for medium values of p\. 
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This fluctuation-induced self- organization is, in fact, observed (see Fig. 2a). The 
entities in subpopulation 2 can agglomerate at sites, where the fluctuations have 
temporarily reduced the density in subpopulation 1 due to disturbances of its 
homogeneous distribution. Later on, subpopulation 1 develops a slightly higher 
concentration at sites where subpopulation 2 clusters. In a similar way, we can 
have drift-induced self- organization for V^ = V > and V 2 = (see Fig. 2b). 
For Di = D 2 = 0, the instability condition (6) is then replaced by 

py 2 (Pi2P2i - Pni^A^Pn + P^ 22 ) 2 > p?p°PnP 22 (n° - V 2 ) 2 , (8) 

where k represents the wave number. (Note that the wave length, which is 
inversely proportional to the respective wave number, is restricted to the values 
A = iAx in our cellular automaton simulations.) 
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Figure 2: Order parameter for the degree of agglomeration as a function of (a) the 
fluctuation strength p 2 = O.OOlpi (with V^ = V 2 = 0) and (b) the drift velocity 
V = V± (with V 2 = = pi = p 2 ) for an example with asymmetric interactions. 
The curves are averages over 20 runs after a time period of 20 000(iVi + iV^) 
update steps for I = 20 cells and Ni = N 2 = 200 entities in each subpopulation. 

For symmetric cases with P n = P 22 = P^0, P\ 2 = P 2 i = Q, and = p 2 = 
p ^ 0, we find Q 2 — P 2 > (V® — V 2 ) 2 /(2kp) 2 , i.e. a finite drift will usually increase 
the threshold for pattern formation. This is different from the effect of diffusion. 
For the symmetric case with V^ = V 2 = and D\ = D 2 = D, the instability 
condition (6) reads p 2 Q 2 > (pP - D) 2 . That is, we expect a "maximum degree 
of self-organization" for D = max(/xP, 0), and a more or less symmetric behavior 
around this point. What does this actually mean? Fig. 3 suggests that this 
statement applies to the order in the system. That is, for increasing fluctuations 
strength we find noise-induced ordering up to D = max(pP, 0), while we have 
disorder at significantly higher fluctuations strengths. 
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Figure 3: As Fig. 2, but for V® = = and symmetric interactions (see the 
black circles in Fig. la for the location in the parameter space), (a) In cases 
of segregation, the order stays about constant for small noise amplitudes, but 
it drops significantly for large ones. In cases of (b) repulsive or (c) attractive 
agglomeration, a suitable noise amplitude pi = p 2 can increase the order in 
the system by more than a factor of 15. The maximum order is reached at 
Pi = Vi ~ max(3P, 0), corresponding to D\ = D 2 = max(pP, 0). Close to 
this maximum, the curves are roughly symmetric, but for less than twice this 
value, the order breaks down completely, (d) The time- dependence of the order 
parameter visualizes the increase of the order in the system. The plot shows 20 
runs for each displayed fluctuation strength, which significantly influences the 
dynamics of the ordering process. After large enough times, we find a typical, 
noise-dependent length scale and level of order in the system. 
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One may think that this fluctuation-induced ordering is due to noise-induced 
transitions from a metastable state (local optimum) to a more stable state of 
higher order (possibly the global optimum). We check this hypothesis for the case 
of repulsive agglomeration, where a coarse-graining appears particularly difficult 
because of the repulsion effect. Our observations are as follows: (i) At moderate, 
but sufficiently large fluctuation strengths, we sometimes observe a "step-wise" 
fusion of agglomerations, which is associated with exponential-like relaxation pro- 
cesses of the time-dependent order parameter to a higher level, (ii) However, the 
main mechanism seems to be that, from the very beginning, fluctuations further 
the formation of larger agglomerations (i.e. suppress the development of small 
ones), which slows down the ordering process in the early stage. (Sometimes the 
system stays disordered for more than 500 000(^ + ^2) iterations and, suddenly, 
the order increases rapidly to a high level), (iii) A careful choice of the noise 
strength can speed up the time-dependent increase of the order very much, (iv) 
After a given, large enough time period, the system has reached a typical level 
of order, which depends significantly on the fluctuation strength. In conclusion, 
one may influence the typical length scale in the system by variation of the noise 
level. 

A variation of the "applied" drift velocity (if possible) or of the fluctuation 
strength together with a proper choice of the "treatment times" would allow one 
to control pattern formation in several respects: (i) the speed of ordering, (ii) 
the typical length scale in the system, and (iii) the level of ordering. A time- 
dependent variation of the control parameters should even facilitate to switch 
between supporting and suppressing structure formation, e.g. between demixing 
and homogenization. In the future, these points may, for example, be relevant 
for the production, properties, handling, and transport of heterogeneous materi- 
als, for flow control [20], and efficient separation techniques for different kinds of 
particles. Due to the general, game-theoretical formulation of the above cellular 
automaton model, its properties are reminiscent of phenomena in physical, bio- 
logical and socio-economic systems. For example, we mention phenomena such 
as the formation of pedestrian lanes (segregation) [17, 18], of ghettos in cities 
(repulsive agglomeration), or of settlements (attractive agglomeration) [21]. 
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